Investigation of target in C-ADS and IAEA ADS benchmark * 

Guojun Hu a - b '\ Hongli Wu a b , Tian Jing a d , 
Xiangqi Wang a,b ' 2 , Jingyu Tang c,a 

"National Synchrotron Radiation Laboratory, USTC 
b School of Nuclear Science and Technology, USTC 
' Institute of High Energy Physics Chinese Academy of Sciences 
''Department of Modern Physics, USTC 



(N 

o 

> 
O 



Abstract 

The spatial and energy distribution of spallation neutrons have an effect on the performance of Accelerator Driven Subcritical 
systems. In this work, the spatial, energy distribution of spallation neutrons and the effect of these factors on proton efficiency 
was studied. When the radius of spallation region increases, backward neutrons were found to have a rather big ratio and have a 
positive effect on proton efficiency. By making better use of these neutrons, we may increase the radius of target to satisfy some 
other requirements in the design of subcritical core. 
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1. Introduction 

China ADS (C-ADS) consists of a high-intensity proton ac- 
celerator in CW (continuous wave) mode with proton energy of 
0.6 GeV in phase three and 1.5 GeV in phase four, a subcritical 
core, and a spallation target of liquid Lead-Bismuth Eutectic 
(LBE or Pb-Bi), which is favored in projects for its low melt- 
ing point, high boiling point and low neutron absorption cross 
section [1]. 

The target should have such a size that it incepts the main 
part of the high-energy cascade and then we get a optimized 
neutron multiplicity [4]. Neutron multiplicity is defined to be 
the number of neutrons produced per beam particle. Increasing 
the size of target will affect the spatial and energy distribution 
of spallation neutrons and thus affecting the performance of the 
total system. So in this work, we discussed the effect of the 
spatial and energy distribution of spallation neutrons. 

Neutron multiplicity, spatial and energy distribution were 
simulated with FLUKA [5, 6] in this work. In these simulations, 
Intra-nuclear cascade, pre-equilibrium, evaporation and fission 
models were activated. Proton efficiency of the system was 
partly simulated with RMC(Reactor Monte Carlo code)[7, 8]. 
IAEA ADS benchmark [9] was an ADS model and was used in 
these simulations, shown in Fig. 2. 

2. Validation of FLUKA 

Double-differential distribution of spallation neutrons were 
measured in the Laboratoire National Sturne experiment [10]. 
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In order to observe the availability of FLUKA in simulat- 
ing spallation reaction, this experiment was simulated with 
FLUKA. Parameters of target, beam and detectors in this sim- 
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Fig. 1 Double-differential neutron production cross-section. 
Data of different polar angle had been scaled with different ratio 
marked in the figure, lcr < 5% 

ulation are listed in Table 1 . Detectors are in a sphere surface 
with radius of R = 800 cm. Detectors are set in polar angle 
of = 0°,10 o ,55 o ,85 o ,130 o ,160°. Double-differential cross- 
section of neutron production is obtained with Eq.(l): 
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Where, dS , dQ. = dS/R 2 , N and L s is area of detector, 
solid angle, atomic density and thickness of target in Table 1. 
dN ns /(N pdm:ir ydSdE) is scored with USRBDX card. 
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Table 1 Parameters of Saturne experiment^ 10] 



Beam 


Target 


Detectors 


Projectile Proton 


Material 


Lead 


R 


800 cm 


E p 1.6 GeV 


As 


1.5 cm 


dS 


1963.5 cm 2 


Direction (0,0,1) 


N 


3.30 x 10 28 nr 3 


dQ, 


3.068 x 10~ 3 sr 


R h 1.4 cm 


L s 


2.0 cm 







Comparison between this simulation and the experiment is 
shown in Fig. 1. Data of the experiment was obtained from 
EXFOR/CRISRS. The result of FLUKA agrees well with the 
experiment. There is difference in high energy due to the defect 
of physical models and the statistic error of simulation. F1UKA 
is suitable to simulate the spallation reaction. 



3. Description of ADS benchmark and proton efficiency 

3.1. Description of IAEA ADS benchmark 

IAEA ADS benchmark was chosen to simulate the effect of 
spallation neutrons. Geometry and parameters of IAEA ADS 
benchmark is show in Fig. 2 and Table 2. In regions "corel" 
and "core2", the ratio between 233-U and 232-Th is 1:9 and 
effective multiplication factor is: k ef{ = 0.96455 + 0.00062. 




Fig. 2 RZ view of IAEA ADS benchmark. The different re- 
gions are:(l)Target; (2)Buffer; (3)Corel; (4)Core2; (5)Core3; 
(6)Reflector; (7)Beam pipe; (8),(10)Moderator 

Liquid Lead is the material of moderator and target in ADS 
benchmark. Spallation region is bounded by the subcritical core 
and beam pipe. We divided the spallation region into "target" 
and "buffer" as used in "Energy Amplifier"[2]. R s and L s is the 
size of spallation region, see Fig. 2. In simulations of FLUKA, 
except for region "target", all other regions were filled with vac- 
uum. 

3.2. Neutron source efficiency and proton efficiency 

Neutron flux <p s in the subcritical core is the solution to the 
inhomogeneous steady-state neutron transport equation Eq.(2): 



Where F is fission operator, A is net neutron loss operator and 
S n is external neutron source. 

The subcritical multiplication factor k s [ 1 1] is defined as Eq.(3): 



< Fcf> s > 



< F(p s > + < S B > 



(3) 



k s describes the number of fission neutrons produced per lost 
neutron in the subcritical system. 

Neutron source efficiency, usually denoted as tp* [11], repre- 
sents the efficiency of the external source neutrons and can be 
expressed as Eq. (4): 
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< F(f> s > is total number of neutrons produced by fission. 

< S n > is total number of external source neutrons. For a given 
values of fc e ff and < Sn >, the larger tp*, the bigger the fission 
power. Definition of neutron source efficiency needs the defini- 
tion of source neutrons < S n >. In this work, < S n > means 
neutrons escaping the spallation region. 

Relationship between kk e s and k s [11] is: 



(i _ _L) = ^(i _ I) 

«keff K s 

Proton efficiency [12], denoted as if/*, is defined to be: 

<Fcf> s > 



<S p > keff-l <S P > 



(5) 



(6) 



A4, S = F</> S +S n 



(2) 



S p is total number of source protons. 

To study proton efficiency of the system, simulations had 
been performed in two steps: 

• protons impact on the target, neutron multiplicity and neu- 
trons escaping the target are recorded with FLUKA. 

• the source neutrons recorded in the first step are then pro- 
cessed as an external neutron source and transported in the 
system. 

In order to reduce the time of simulation, the height and diam- 
eter of the total system was cut to be 360 cm. 



4. Spatial and energy distribution of spallation neutrons 

4.1. Description of spatial and energy distribution of spallation 
neutrons 

Spatial and energy distribution of spallation neutrons have an 
effect on proton efficiency. 232-Th,238-U, Am can be used as 
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Table 2 Parameters of ADS benchmark[9] 

The element of IAEA ADS benchmark at BOL (20 °C) 10 24 crrT a 



Nuclides 


Corel 


Core2 


Core3 


Reflector 


Moderator, Target, Buffer 


232-Th 






7.45 






232-Th,233-U 


6.350 


7.450 








O 


12.70 


14.90 


14.90 






Fe 


8.100 


8.870 


8.870 


6.630 




Cr 


1.120 


1.060 


1.060 


0.800 




Mn 


0.046 


0.051 


0.051 


0.038 




W 


0.046 


0.051 


0.051 


0.038 




Pb 


17.70 


15.60 


15.60 


0.024 


30.50 



The source of IAEA ADS benchmark 


Projectile 


Proton 


Kinetic Energy:£p 


1.0 GeV 


Beam radius:/?b 


10.0 cm 



fuels of ADS. An external neutron hardly ever induce fissions 
of these fuels when its kinetic energy is below 1 MeV[12]. The 
energy spectrum of neutrons in subcritical core is different from 
that of spallation neutrons; the former is determined by geom- 
etry and composition of the system, while the later is mostly 
determined by the size of spallation region. The spatial distri- 
bution of spallation neutrons determines the distribution of neu- 
tron flux in the subcritical core and thus determining the proton 
efficiency of the system. In this work, spatial and energy distri- 
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Fig. 3 RZ view of target model used in FLUKA. The dif- 
ferent regions are: (l)Target; (2)Subcritical core; (3)Beam 
pipe; (4)Moderator; (5)Void. In simulations of FLUKA, region 
Moderator and Subcritical core is filled with vacuum. Spal- 
lation neutrons include : Backward(l— »3,4); Forward(l— »5); 
Radial(1^2); Survival(4^2); Leak(3,4^5) 

bution of spallation neutrons was simulated with FLUKA. The 
target model used in FLUKA is shown in Fig. 3. Neutrons are 
divided into forward, radial, backward neutrons according to 



the surface the neutrons cross when they escape the target. Ra- 
dial neutrons enter the subcritical core directly. Forward neu- 
trons hardly ever enter the subcritical core. Some backward 
neutrons make it to the subcritical core through the moderator 
between the beam pipe and the subcritical core (region "10" in 
Fig. 2), called survival neutrons; other backward neutrons enter 
the reflector or leak out from the beam pipe, called leak neu- 
trons. A clear explanation of these neutrons is shown in Fig. 3. 
Except for region "target", all other regions are filled with vac- 
uum. 

As for spatial distribution, what we concern here was the neu- 
tron current surface density / r (Z) in the side surface of target. 
Z is the distance between a point in the side surface and the 
top surface of target. We also cared about the number of back- 
ward, forward and radial neutrons. As for energy distribution, 
the energy spectrum of radial neutrons were mostly discussed. 

4.2. Spatial and energy distribution of radial neutrons in C- 
ADS 

Radial neutrons are the main part of spallation neutrons and 
their spatial distribution is closely related to the neutron flux 
in the subcritical core. So we studied this factor quantitatively 
with F1UKA. Remember that we mainly care about neutron cur- 
rent surface density here. 

In this simulation, proton source is a circular uniform beam 
with radius of R b = 15 cm in C-ADS HEBT (High Energy 
Beam Transport ) [13]. Simulation had been performed with 
parameters listed in Table 3. To obtain the neutron current sur- 
face density, subroutine mgdraw.f was modified and activated 
to record neutrons escaping the "target" region in Fig. 3. Then 
we constructed the spatial and energy distribution of neutrons 
with the information recorded. The result of this simulation is 
shown in Fig. 4. Average kinetic energy of spallation neutrons 
in different position was used in calculating energy distribution: 
E r (Z). As can be seen in Fig. 4, as Z increase: E r (Z) reaches a 
minimum and then saturates; 7 r (Z) reaches a maximum and then 
decreases to zero. Based on this result, if the position of max- 
imum 7 r (Z) locates at the axial center of subcritical core(that 
is z — cm), neutron flux <f>$(z) in the subcritical core will be 
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Table 3 Parameters of target and beam of C-ADS [13] 



Target characteristic 


Beam 


Material 50%Pb,50%Bi 


Projectile Proton 


R s 25.0 cm 


E p 1.5 GeV 


L s 100.0 cm 


Source HEBT 



symmetric (with axis at z — cm). Generally, an symmetric dis- 
tribution of 4> s (z) is a key factor in designing the system. In this 
way, we may optimize the position(relative to the subcritical 
core) of spallation target, discussed in Sec. 4.4.1. In this case, 



Spatial and energy distribution of radial spallation neutrons 
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Fig. 4 Neutron current surface density / r (Z) and average kinetic 
energyE^Z) as a function of Z. / r (Z) is normalized to per source 
proton. 

E r (Z) changes slightly, (E (rmax ) - £rmin)/£(raverage) ~ 0.141. En- 
ergy distribution of radial spallation neutrons has a influence on 
the neutron flux of the subcritical core, but it is a less important 
factor than the spatial distribution of radial neutrons. Detailed 
effect of the energy distribution needs further research. 

4.3. Backward neutrons and radial neutrons as a function ofR s 

The size of subcritical core changes with the radius of spalla- 
tion region. It was proved that proton efficiency decrease when 
radius of spallation region (called target in [12]) increases due 
to the soft and leakage of neutrons [12]. Changing the size of 
the subcritical core will bring to change of several factors that 
influence the proton efficiency. In this work, as is shown in 
Fig. 5, a "void" region is inserted between the spallation region 
and the subcritical core to study the effect of backward neu- 
trons. Changing the thickness of this "void" region, we could 
simulate different R s without changing other parameters, such 
as the size of the subcritical core. This "void" region does not 
exist in practical design of ADS, but it is useful to study the 
effect of backward neutrons. 

To study the effect of backward and radial neutrons on proton 
efficiency with ADS model, parameters of IAEA ADS bench- 
mark were used, as listed in Table 2; L s = 100 cm. We 
use the target model in Fig. 3 to study the number of back- 
ward neutrons and the soft of radial neutrons when R s changes. 
Four different thickness of "Void" region were chosen to get 
R< = 12, 18, 24, 30 cm. The result is listed in Table 4. Remem- 




Z/cm 



Fig. 5 RZ view of modified ADS benchmark. The different 
regions are:(l)Target; (2)Buffer; (3)Corel; (4)Core2; (5)Core3; 
(6)Reflector; (7)Beam pipe; (8)Moderator; (9) Void. 

ber that in all simulations with FLUKA, region subcritical core 
is filled with vacuum. 

Table 4 Characteristic of spallation neutrons, lcr < 1.0% 
Number of spallation neutrons in different direction 



R s /cm 


12 


18 


24 


30 


Radial 


20.665 


21.838 


21.219 


20.199 


Forward 


0.007 


0.034 


0.110 


0.263 


Backward: 


3.729 


6.368 


8.702 


10.767 


Survival 


3.530 


5.697 


7.270 


8.281 


Leak 


0.199 


0.671 


1.432 


2.486 


Ratio of high-energy Radial neutrons 


> 100MeV% 


1.70 


0.99 


0.62 


0.41 


> 20 MeV% 


1.26 


0.86 


0.58 


0.39 


> 1 MeV% 


59.62 


48.3 


38.97 


31.04 



In Table 4, as R s , we find that: (1) the number of radial neu- 
trons changes slightly and the spectrum of radial neutrons be- 
comes softer (ratio of neutrons with kinetic energy higher than 
1 MeV decrease from 59.62% to 31.04%); (2) the number of 
backward neutrons increases greatly (from 3.729 to 10.767); (3) 
the number of survival neutrons increases considerably (from 
3.530 to 8.821). Explanation to these results may be: most neu- 
trons, especially high energy neutrons are produced near the 
top surface of target; as R s increases, it is more likely that these 
neutrons make it to the top surface of target and then leak out 
in backward direction. In conclusion, when the radius of spal- 
lation region increases, number of survival neutrons becomes 
considerably big, which will contribute a positive effect to pro- 
ton efficiency; radial neutrons become softer, which will con- 
tribute a negative effect to proton efficiency. Comparison of 
these two effects are discussed below. 

In this simulation, material of region "moderator" in Fig. 3 
was vacuum, not Lead; replacing the vacuum with Lead would 
bring a slight decrease to the number of survival neutrons 
due to the scattering and absorption of neutrons in the "mod- 
erator". Another simulation had been performed; for R s = 
12, 18,24,30 cm, survival neutrons decrease to 3.409, 5.098, 
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6.162, 6.838 per source proton. This difference does not 
severely affect the fact that a rather big number of backward 
neutrons will enter the subcritical core when the radius of spal- 
lation region is big. 

4.4. Effect of spatial and energy distribution of spallation neu- 
trons 

4.4.1. Optimizing target position with spatial distribution of ra- 
dial neutrons 

The spatial distribution of radial neutrons may be used to op- 
timize the axial position of target, as discussed in Sec. 4.2 . In 
this work, target position is determined by the top surface of 
target Zt op . Optimization of Z top had been performed with the 
following considerations: 

• L s will be larger than the range of proton in Lead, about 
100 cm for 1.0-1.5 GeV protons. In this way, I r (Z) and 
neutron multiplicity will not change when Z top changes. 

• The bottom of target, Zbottom, remains at 75 cm; Z top 
changes from -60 cm to -15 cm and L s varies from 90 cm 
to 135 cm. 

• All other parameters do not change, such as R s = 32.5 cm 
and composition of the subcritical core. k e ff has a little 
change and will be considered in the calculation of proton 
efficiency. 

• L s reaches its saturation, so< S n > /< S p > does not 
change when Z top changes. <p* is equal to if/*. 

The position of maximum 7 r (Z) varies with the radius of spal- 
lation region and the energy of protons. For ADS benchmark, 
Rs = 32.5 cm, L s = 100 cm, E p = 1.0 GeV, the maximum 
reaches at Z a 20 cm. So we expected a maximum proton ef- 
ficiency when Z( t0 p) = -20 cm. Parameters of this simulation 
are listed in Table 3. cp* and s (z) was obtained as Z top changed. 
The result is shown in Fig. 6 and Fig. 7. 
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Fig. 7 Neutron flux <p s (z) as a function of Z, 



Ion 



In Fig. 6, we find that both £ s and tp* reaches a maximum 
value at Z top = -25 cm. This result is a bit smaller than the po- 
sition, -20 cm, we expected. This difference may be explained 
by the fact that the distribution of / r (Z) is not perfectly symmet- 
ric, there is a "tail" at high Z. In Fig. 7, as Z top departs from 
-25 cm, the symmetry axis of (f> s (z) departs from z — cm. 
This deviation will increase the leakage of neutrons from the 
subcritical core to the reflector and thus decreasing the proton 
efficiency. 

Survival neutrons also influence <p*. The number of survival 
neutrons will decrease when Z top becomes more negative and 
this will bring a decrease to tp*. The effect of 7 r (Z) and survival 
neutrons are coupled. But when Z top becomes more positive, the 
effect of survival neutrons becomes smaller. In conclusion, both 
(f> s (z) and tp* support the assumption about target position. The 
effect of survival neutrons makes the assumption less persuasive 
and will be discussed below. 



4.4.2. Effect of radial and backward neutrons on proton effi- 
ciency 

Change of R s will bring a change to the characteristic of ra- 
dial and backward neutrons as discussed in Sec. 4.3. These two 
factors are coupled. Based on simulations in Sec. 4.3, proton 
efficiency if/* of the modified ADS benchmark was simulated 
for different R s . Parameters of these simulations are listed in 
Table 2, the result is listed in Table 5. Changing the thickness 
of void region brought to a slight change of k e s- This change 
was considered in the calculation of if/* . 

Table 5 Proton efficiency varies with /? s , Icr < 1.0% 



Proton efficiency R $ 


R s /cm 


12 18 24 


30 


k eS + std 


0.94828 0.94886 0.94972 


0.95148 


0.00026 0.00026 0.00025 


0.00026 


r 


26.803 27.499 27.882 


28.002 



Fig. 6 Neutron source efficiency tp*, subcritical multiplication 
factor k s and effective multiplication factor k e s as a function of 

Ztop 



We find that: as R s increases, if/* increases slightly. This re- 
sult seems paradox to the conclusion that the bigger of target 
radius, the smaller of proton efficiency in [12]. In fact,as is 
seen in Table 4, it is the trade-off between the soft of radial 
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neutrons and the increase of survival neutrons. Source effi- 
ciency of backward neutrons in inducing fissions is lower than 
these radial neutrons, but the rather big number of survival neu- 
trons neutralizes the negative effect of radial neutrons. This 
result does indicate that the effect of backward neutrons is pro- 
nounced. When radius of spallation region increase, number of 
neutrons produced per beam particle does increase, but these 
extra neutrons are not effectively used. Based on this simula- 
tion, it is difficult to draw the conclusion that the radius of spal- 
lation region should be larger; but if we need to increase radius 
of spallation region to satisfy other requirements in designing 
the system, such as decreasing the damage of high-energy spal- 
lation neutrons to structural material, we may try to make better 
use of these backward neutrons and avoid severely affecting the 
proton efficiency. At present, we may replace the moderator 
between the beam pipe and the subcritical core with fuels and 
increase the radius of spallation region; or we may design an an- 
nular target (ratio between area of top surface and side surface 
is smaller than that of a cylinder target) to reduce the number of 
backward neutrons. 

5. Conclusion 

In this work, spallation region and subcritical core in ADS 
are combined to study effect of the spatial, energy distribution 
of spallation neutrons. The characteristic of radial and back- 
ward neutrons were mainly discussed. 

As for radial neutrons, neutron current surface density I r (Z) 
was simulated with FLUKA and used to optimize the position 
of target. An assumption about this optimization was made and 
tested. As for the design of C-ADS, the considerations made in 
Sec. 4.4. 1 may be used to optimize position of target. As for 
backward neutrons, survival neutrons were mainly discussed 
and found to have a considerable ratio when radius of spalla- 
tion region increases. These survival neutrons have a positive 
effect on the proton efficiency. If we make better use of these 
survival neutrons, we may increase the radius of spallation re- 
gion to satisfy other requirements in the design of subcritical 
core without affecting proton efficiency severely. 
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